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The present status of the heavy-ion program to search for quark-gluon plasma 
is reviewed. The goal of this program is to recreate the Big Bang in the labo- 
ratory, by generating small chunks of exploding quark-gluon plasma ("The Little 
Bang" ) . I argue that the analogues of the three pillars of Big Bang Theory (Hubble 
, flow, microwave background radiation, and primordial nucleosynthesis) have now 

*^ 1 been firmly established in heavy-ion collisions at SPS energies: there is convincing 

\I ' evidence for strong radial flow, thermal hadron emission, and primordial hadrosyn- 

Cn ' thesis from a color-deconfined initial stage. Direct observation of the quark-gluon 

\| ' plasma phase via its electromagnetic radiation will be possible in planned collider 

Cn experiments at higher energies. 

o' 

5^ \ 1 Introduction: What Are We Looking For? 



f^ • The highest particle and energy densities existed for a fleeting moment in the 

early history of our universe, shortly after the "Big Bang" . In heavy-ion col- 
lisions at ultrarelativistic energies one hopes to be able to recreate such high- 

(-H \ density matter and to study its properties. At presently available beam enen- 

gies the goal is to pass the energy density threshold for color-deconfinementid, 
Ecr 5; 1 GeV/fm'^ . If such energy densities can be reached and the energy 
sufficiently thermalized, strongly interacting matter will manifest itself as a 



H \ quark-gluon plasma: the hadronic constituents (quarks and gluons) become 

deconfined, and hadrons loose their identity. 

1.1 The Big Bang model as an example 

How do we know that high energy densities existed in the Early Universe? 
How is the Big Bang Theory "proved"? It is worthwhile to contemplate this 
question for a minute. According to the textbooks the Big Bang Theory rests 
on three pillars: (1) The Hubble law relating galaxy distances to their reces- 
sion velocities, which proves that our universe is expanding; (2) the cosmic 
microwave background (CMB) with its uniform temperature of 2.7 K; and (3) 
the successful prediction of the primordial abundances of small atomic nuclei 
during nucleosynthesis. The present expansion rate of our universe reflects 



its initial conditions and the equation of state of the matter in the uni- 
verse, integrated over its lifetime. Primordial nucleosynthesis happened when 
the universe had cooled down to about T ~ 100 keV, about 3 minutes after 
the Big Bang; this is the stage of chemical freeze-out of our universe at 
which its original chemical composition was fixed (later to be modified inside 
stars by non-cosmological mechanisms). The microwave background radiation 
signals the recombination of atoms and the depletion of free charges in the 
early universe, rendering it transparent to photons. This thermal freeze-out 
of our universe happened much later, about 400000 years after the Big Bang 
when it had cooled down to T ~ i eV ~ 3000 K. 

Before thermal freeze-out, the universe was in nearly perfect local and 
global thermal equilibrium (except for the neutrinos which froze out earlier 
and today have a temperature of 1.95 K). Before the onset of nucleosynthesis, 
it was also in a state of chemical equilibrium. Thermodynamic equilibrium 
being a state of maximum entropy or minimal information, all memory about 
the earlier stages of the universe got wiped out, leaving practically no directly 
observable traces. In particular, the color-confining phase transition with its 
formation of hadrons from quarks and gluons at about 10 /iS after the Big 
Bang, which we now want to study in the laboratory, is in cosmology safely 
hidden behind the curtain of thermal CMB photons. No direct probes from 
the early stages are known which survived this tendency for re-equilibration; 
given the roughly 18 orders of magnitude between the (slow) expansion rate of 
the universe and the (fast) strong, weak and electromagnetic interaction rates 
among the particles this is, of course, not surprising. At the moment our best 
bet for direct signals from the very early universe seems to be a detailed and 
quantitative analysis of the fluctuations in the CMB spectrum which occur 
at the level of about 10~^ of the average signal. 

Still, with only a single observed event and little or no direct experimental 
evidence about its very early stage to show, the Big Bang Theory has become 
the almost unanimously accepted model for the evolution of our cosmos. How 
can this be explained? The most important reason is, I think, the fact that 
we can analyze our present state (the final state of the "Big Bang experi- 
ment") by making exceedingly accurate observations. This provides stringent 
constraints for the extrapolation backwards in time. Furthermore, general rel- 
ativity, which has been tested elsewhere, provides a robust and highly reliable 
theory for the dynamical evolution of the cosmos, given the equation of state 
of the matter. And the latter again is known experimentally and theoretically 
to be simple during most parts of the dynamical evolution, with the exception 
of a few critical stages (typically near phase transitions) which correspondingly 
draw most of the present theoretical attention. 



1.2 What is a quark- gluon plasma? 

The initiated reader will not have missed the strong parallelism between the 
above discussion of the Big Bang and the field of relativistic heavy-ion colli- 
sions. Very similar concepts (indicated above by emphasized key-words) play a 
crucial role on both contexts. Why, then, do we not have a generally accepted 
theory for the "Little Bang"? 

The main reason is the small geometric size of the object under study. 
Contrary to the Early Universe, the region of high density is limited to a re- 
gion of at most a few 10 fm in diameter, surrounded by vacuum. Once the 
energy in the reaction zone begins to thermalize, it builds up pressure which 
leads to explosive expansion of the fireball. The associated expansion time 
scale is no longer 18 orders of magnitude, but only perhaps a factor 10 to 
100 larger than that of the microscopic equilibration processes. This intro- 
duces basic uncertainties into dynamical evolution models which to eliminate 
requires hard work and many detailed tests. It is fair to say that the goal of a 
generally accepted "standard model of nuclear fireball dynamics" has not yet 
been reached. 

Things would be much simpler if, as in cosmology, we were sure that 
hydrodynamics can be used. In fact, the conditions for creating a quark-gluon 
plasma (QGP) and for being able to describe its dynamics hydrodynamically 
are very closely related. Both involve an assessment of the degree of local 
equilibrium that can be reached. I would list the following criteria for being 
permitted to call a small, short-lived, rapidly expanding system of quarks and 
gluons a quark-gluon plasma: 

1. Sufficient energy density for deconfinement, e 3> Ccr — 1 GeV/fm^. 

2. Locally thermal momentum distributions with the same temperature T 
for all particle species, in a limited momentum range 0<p < 6 T. To 
establish collective plasma properties and to validate a hydrodynamic 
approach, it is not necessary to fully equilibrate the tails of the momen- 
tum distributions. In fact, given the short lifetime and small size of the 
reaction zone, some fast particles produced by hard QCD processes will 
always escape, giving a power-law tail to the momentum spectrum. Note 
that chemical equilibrium, is also not necessary: the relative abundances 
of gluons and light and heavy (anti)quarks need not be equilibrated for 
plasma properties and hydrodynamic behavior to manifest themselves. 
Chemical non-equilibrium requires, however, the introduction of addi- 
tional chemical potentials in the equation of state. 



3. SufRcicntly large volume and strong rescattering. This and the foUowing 
criterium are the most subtle ones: the relevant volume is characterized 
by the mean free path which again reflects the scattering rate. To ensure 
thermalization by kinetic equilibration we must require 

If » 1 (1) 

where in perturbation theory the mean free path is givenl3l3 in terms of 
the time-dependent temperature T{t) ^ e(r)^/'* bypl 

Amfp(T) = rscatt(T) (v) ^ 2f \: n / t ^^ (^) 

#T(r)a^(r)ln(l/a(r)) 

(with # indicating a number proportional to the number of degrees of 
freedom in the plasma), and Vhom(p) is the "homogeneity volume" of 
particles with momentum px 

Vhom{m±, y,T) = 2tt Ri\{m±,y,T) R']_{m±,y,T). (3) 

-R|| , R±_ are the longitudinal and transverse "homogeneity lengths" as 
measured by Bose-Einstein correlations. For thermalized exnanding fire- 
balls they can be estimated by simple analytical expressionsa which, as- 
suming boost-invariant longitudinal expansion^ lead to 

TT / \ V'goom(''"j 

Vhom[m±,y,T) = 






— « ^geom(r) (^^ j . (4) 

Here the limit is for large m±^ , rjf (r) is the average transverse flow rapidity 
of the fireball, A?7(t) its longitudinal extension in space-time rapidity, 
and Vgcoi-a{T) = 7ri?^(r) • 2tA77(t) is the geometric fireball size (growing 
with time). The interesting ratio (|l|) then becomes 



.,, , - ^gcom {t)T\t) ■ #^a^(r) ln^(l/a(r)) [^^\ . (5) 

Entropy conservation requires the first factor to be time-independent 
whereas the second factor grows with time. One sees that the cri- 
terium (Q) is the more easy to fulfill the higher the initial temperature 



"This implies a transport cross section iTtrans ~ a^ \n(l/a)/T^ which grows as T decreases. 



(To = r(To)), the larger the number of massless degrees of freedom m the 
plasma (#), the more strongly they interact (a), and the smaller their 
transverse mass {m±). Particles with large m± will always be difficult 
to thermalize. 

4. Sufficiently long lifetime: 

Tscatt ^ -^^^ < 6*"^ = = Tcxp. (6) 

c o ■ u 

The divergence of the flow velocity field 9{x) — d ■ u{x) gives the local 
expansion rate in the fireball; it is the analog of the Hubble constant 
in cosmology. For d-dimensional scaling expansion it is given by 0~^ = 
Tcxp = rZ^whereas the temperature drops like T{t) = ro(ro/r)^/'^, with 
TqTq w iH'EI Putting this together we find 



Tscatt dinrY'^-^ 



#a2(r)ln(l/a(r))- 



(7) 



For d—\ the expansion looses and the conditions for thermalization get 
better and better with time; the case d—2> is marginalB, and only the 
increasing coupling strength cxLt) at lower temperatures works in fa- 
vor of increased thermalizationB Of course, after hadronization the cross 
sections become r-independcnt while the particle density continues to 
decrease such that eventually thermalization stops and the particle mo- 
menta freeze out. 

Clearly, the last two criteria require sufficiently strong coupling a. Small values 
of a are certainly disadvantageous for thermalization in heavy-ion collisions; 
fortunately, in real life a seems to be sufficiently large (much to the dismay of 
the practitioners of perturbative QCD). 

2 Heavy-Ion Observables: Hard and Soft Probes 

The observables in relativistic heavy-ion collisions can be divided into two 
classes: hard and soft probes. Hard probes are created early in the collision 
and, due to the finite size of the reaction fireball and their relatively small rein- 
teraction cross section, they decouple early. They include hard direct photons 
and lepton pairs, hadronic jets, and charmonia {J/ipyip', . . .). Soft probes are 
the light quark flavors u, d, s and the corresponding hadrons, as well as the soft 
electromagnetic radiation emitted by them. They are created throughout the 
collision history, and the strongly interacting ones decouple late, triggered by 



the expansion of the fireball rather than by its finite size. Soft probes include 
the hadron abundance ratios, their niomentum spectra and their momentum 
correlations. 

The soft probes contain equivalent information as the known signatures 
of the cosmological Big Bang: Hadron abundancies give information about 
the chemical equilibration time scales; here strangeness plays a crucial role 
since the time scale t^s for creating strange quark pairs is roughly of the same 
order of magnitude as the fireball lifetime such that moderate changes in t^ 
due to new physics (e.g. color deconfinement and the onset of gluon fusiorJH) 
can cause major effects. The hadronic momentum spectra and two-particle 
correlations provide information on thermalization and collective flow. Since 
flow is generated by pressure, the observed flow pattern in the final state gives 
a time integral of the pressure history in the fireball (and thus its equation 
of state), with different types of flow (radial expansion in central collisions, 
directed and elliptic flow in semipedpheral ones) receiving different weights 
from early and late dynamical stagesil The observed flow patterns are thus the 
"Little Bang" analogues of the Hubble expansion in cosmology. 

The Little Bang can therefore be reconstructed from the soft probes in very 
much the same way as the cosmological Big Bang; hard probes, which escape 
directly from the early stages of the Little Bang and for which no cosmological 
analogue exists (due to the absence of spatial boundaries of the universe), can 
then be used to check the consistency of the reconstruction. The presently 
available data allow for the realization of the first part of this program; the 
second part will require the higher collision energies becoming available at 
RHIC and LHC in the next few years. 

3 Hard Probes: The Present Status 

An extensive review of the present status of QGP signatures has just appeareda 
ft contains a comprehensive list of original references to which f refer the reader 
for details. 

Direct photons of other than hadronic decay origin have been searched for 
by the WA80, WA98 and CERES/NA45 coUaborations at the CERN SPS in 
S+Au and Pb+Pb collisions. So far only an upper limit of about 5-7% above 
hadronic decay background could be established in S-f Au collisions; the analy- 
sis of Pb-|-Pb data is still in progress. This eliminates hydrodynamic expansion 
models with equations of state which do not include a phase transition to QGP 
or at least a strong softening of the equation of state by resonance and string 
excitations at high temperatures. A positive signal of thermal QGP radiation 
(whose rate grows with T*) is, however, only expected at higher beam energies 
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where larger initial temperatures can be reached. Whether it can be extracted 
from the hadronic backgimmd depends on the unknown transverse momentum 
distribution of the latterEHl (which is affected by collective expansion). 

An excess of dimuons in the mass region below the J /ip has been reported, 
in one form or another, by the HELI0S3, NA38 and NA50 collaborations in S- 
and Pb-induced heavy-ion collisions. This excess appears to have a non-linear 
dependence on the charged hadron multiplicity (collision centrality) . Its origin 
is presently unclearli^ 

The CERES/NA45 collaboration has seen in S-f Au and Pb-f-Au collisions 
a strong excess of e+e" pairs with masses 200 MeV< rriee < rnp, with a nop^ 
linear multiplicity dependence and concentrated at low transverse momentumlla 
It can probably be explained by hadronic mechanisms, but requires a very 
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Figure 1: CERES dilepton mass spectralid and theef^tical simulations with an expanding 
fireball model, using different theoretical approachesE3 to the in-medium dilepton rates. 



dense and hot hadronic medium in which pions and p-mesons violently rescat- 
ter, leading to a broadening of the vector meson spectral densitieslillFig. |l|). 



Jet production and jet quenching by the dense fireball mediumEzI cannot 



be studied at present SPS energies due to insufficient rates, but will become 
accessible at RHIC. 

Much recent attention focussed on the discovery by NA38/NA50 of "ano- 
malous" suppression of J/^ and ijj' vector mesons in central nucleus-nucleus 
collisions. Different collision systems can be compared by introducing a theo- 



retical auxiliary variable, the average path-length L which the cc-pair must 
traverse before escaping the reaction zoneEj "Normal" suppression is defined 
by an exponential attenuation Ba^/aoY oc exp(— pCTabs^) (straight line in 
Fig. y) with the normal nuclear density p and a fitted ipN absorption cross 
section (Tabs- A deviation from this behaviour is first observed for the "0' near 
L = 5 fm, then for the J/ip near L = 7.5 fm. Since about 32-40% (5-8%) of the 
observed J/-0 stem from radiative Xc {4'') decays, it is not unlikely that the 
drop near i = 7.5 fm indicates anomalous Xc suppression, and that anomalous 
suppression of the J /^ itself requires even larger values of L. It is interesting 
that this supression pattern follows the binding energies of the corresponding 
charmonium states: suppression of the more strongly bound states requires 
a higher density and/or lifetime of the fireball, here parametrized (perhaps 
not very fortunately) via L. At present the suppression mechanism is not 



— s 






i 7« 


- 


K U 


O pt450ae\7cl-A iA^p.d) 


y, w 


~^jl^ « p(200GcV/ci-A (A = W,Ui 




•f^^-^ ■ .S' (32 x200GeVk)-U 


^ 40 


- ^^^^ 1 • ['b (2lifi K ISX GeV/i:j - Pb 


"^ 30 


^^"^^1 










=0 


rf--. 


M 


' --, 




By.yirsU/'Vi/cfDY) ^Ae-fi^-bf- ^j!^A 




a.i,, =S.8+/-0.6mb H 




p„ =0.17fm-' 


10 


-| 1 1 1 




2.5 5 7.5 10 


L(fm) 



p (200 GeV/c) - A (A=p,tl) 
p (450 GeV/c) ■ A (A = W,U) 
S (32 X 200 GeV/c) - V 
Pb (208 X 158 GcV/c) - Pb 




Figure 2: Preliminary data on {J/i))/DY and i\)' jDY from the NA38/NA50 collaborationlla, 
plotted as a function of the Glauber model nuclear thickness parameter L (see text). 

fully understood theoreticallytj but one condition appears to be unavoidable: 
strong rescattering of the cc-pair in a very dense environment, probably of 
partonic origin. A consistent quantum mechanical and dynamical theory of 
this phenomenon is urgently needed. 



4 Soft Probes: Hadronization, Thermalization and Flow 

Soft hadrons, which constitute the bulk of the produced particles, provide a 
much richer body of experimental information which by now has been ana- 
lyzed in considerable quantitative detail. As I will show, they give convincing 
evidence that we have seen "the Little Bang" . 



4-.1 Primordial hadrosynthesis 



It has been known for many years that hadron production in high energy 
physics exhibists striking statistical features. Recently several rather detailed 
analyses were performed of the relative yields of different hadronic species 
produced in e'^e^, pp, pp and AA collisions, with the question in mind pa 
what extent the hadronic final state reflects a state of chemical equilibriumxJi 
The result is shown in Fig. S 
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Figure 3: Compilation of freeze-out points in e+e" collisions and in heavy-ion collisions from 
SIS to SPS energies. Filled symbols: chemical freeze-out points from hadron abundances. 
Open symbols: thermal freeze-out points from momentum spectra and two-particle correla- 
tions. The shaded region indicates the parameter range of the expected transition to a QGP. 

The single e+e~ point represents a large number of e'^e^ , pp and pp colli- 
sion systems atiiifferent center of mass energies 20 GeV < -y/s < 900 GeV which 
were all founclE£l to reflect chemical equilibrium hadron abundancies at a uni- 
versal hadronization temperature Thad ~ 175 MeV. The only clear deviation 
from chemical equilibrium in these systems is an undersaturation of overall 
strangeness, reflected in the ratio of produced strange to non-strange quark 
pairs As =: 2(ss)/((uu) + (d(i))| produced ~ 0.2—0.25, again almost independent 
of y^EjO. But even if the total number of ss valence quark pairs is below the 
chemical equilibrium value, the available s and s quarks are distributed among 
the various strange hadron species according to the law of maximum entropy. 

Given the small size of these collision systems it is impossible to imagine 



that this can be the resuh of cquihbration by hadronic rescattcring. It must 
reflect pre-established statistical equilibrium, i.e. the statistical filling of the 
available hadronic phase-space according to the law of Maximum Entropy at 
the point of hadron formation. The temperature Thad should be understood 
as a Lagrange multiplier for the energy density at which this happens; its 
universality and numerical value tell us that hadronization always happens at 
the same critical energy density Ccr ~ 1 GeV/fm'^ . At higher collision energies 
hadrons are not produced, at higher temperatures, but over a larger volume at 
the same energy densit'i)L3. After formation the hadrons do not rescatter, and 
the observed hadron abundances thus reflect the primordial values established 
at the point of hadronization. 

It is interesting that in heavy-ion collisions at the SPS (S-l-S, S-|-Ag, 
Pb+Pb) one finds again chemical equilibrium hadmn abundances at the same 
chemical freeze-out temperature of about 180 MeVta This suggests hadron for- 
mation by the same statistical phase-space filling mechanism as in e'^e~ and 
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Figure 4: The strangeness suppression factor Aa as a function of v^tJ- The two points each 
for pp coUisions reflect the inclusion (exclusion) of the initial valence quarks. 

pp collisions (reflecting the Maximum Entropy Principle), followed by imme- 
diate freeze-out of the hadron yields. Since changes in the relative abundances 
require inelastic collisions whose cross sections are mostly quite small, this is 
not unexpected, especially if the fireball features collective expansion (which 
is known to foster freeze-out) already at hadronization. 

What is different from elementary collisions is that in heavy-ion collisions 
a considerabljj-Jarger global strangeness fraction As ~ 0.4 — 0.45 is mea- 
sured (Fig. H)tj Following the above argument and realizing that hadronic 
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strangeness- creating cross sections are particularly smallpl this can only be 
due to processes before hadronization. Indeed, all known microscopic kinetic 
models based entirely on pp input and hadronic reinteraction dynamics fail to 
reproduce this strangeness enhancemeniz3 Heavy-ion reactions thus generate a 
prehadronic state with different dynamics (presumably a longer lifetime) than 
in e+e~ and pp collisions. The factor 2 increase in the strangeness fraction 
As indicates a shnrt strangeness saturation time scale r^j in this state, as pre- 
dicted for a QGPlJ In fact, SoUfrankEd has argued that the observed strangeness 
enhancement in A + A collisions at the SPS is consistent with the hadroniza- 
tion of a fully equilibrated QGP at the hadronization temperature Thad if the 
hadronization process itself leaves the strangeness/entropy ratio unchanged. 
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Figure 5: Enhancement factor for the mid-rapidity yields per participating nuclncpji in Pb+Pb 
relative to p+Pb coUisions for various strange and non-strange hadron speciesE3 

The global strangeness enhancement reflected in A^ occurs_already in S+S 
collisions and remains roughly unchanged in Pb-t-Eh. collisionsEj A recent com- 
parison of p+Pb and PbH-Pb collisions at the SPScJ (see Fig. ^ further shows 
that while the bulk of the strangeness enhancement is carried by the kaons and 
hyperons (A, E), which are enhanced by about a factor 3 near midrapidity, the 
enhancement is much stronger for the doubly and triply strange baryons S 
and fl and their antiparticles, with an enhancement factor of about 17 (!) for 
n + n at midrapidity. Such a rise of the enhancement with the strangeness con- 
tent of the hadrons is natural in statistical hadronization models, but contra- 



il 



diets expectations based on the respective production thresholds in hadronic 
(re)interactions. WA97 also founcO that all of the enhancement factors in 
Pb+Pb are independent of the collision centrality, i.e. of the size of the midra- 
pidity source, from about 100 to 400 participating nucleons, and that similar 
enhancement patterns are seen in S+S collisions^ Thus, whatever causes the 
enhancement in central Pb+Pb collisions exists already in S+S collisions! 

The existence of a prehadronic stage without color confinement, both 
in S+S and Pb+Pb colhsions at the SPS, is supported by^ recent argu- 
ment by Bialascj which generalizes similar ideas by RafelskiE3 by removing 
the assumption of thermal equilibrium. Bialas points out that by considering 
baryon/antibaryon production ratios the unknown baryon wave function drops 
out, and one can very easily test whether or not the baryons were formed by 
statistical hadronization (coalescence) of uncorrelated quarks. He finds that the 
data from S+S and Pb+Pb, but not those from p+Pb follow the corresponding 
simple "quark counting rules" ; the conspicuous absence of correlations among 
the quarks is intjexpreted in terms of a color-deconfined initial state in the 
nuclear collisions Eil 

While at the SPS chemical freeze-out appears to happen immediately at 
hadronization (such that the measured hadron yields reflect the abundances 
from the primordial hadrosynthesis) , Fig. || suggests that at the AGS and SIS 
chemical freeze-out occurs at lower temperatures, after some further "chemical 
cooking" . This is probably due to longer lifetimes and slower expansion of the 
reaction zone at lower beam energies. 



4-2 Thermal hadron radiation and radial flow 

I now discuss the analogues of Hubble expansion and the cosmic microwave 
background. In heavy-ion collisions the latter is not generated by photons (be- 
cause these decouple immediately after their production, i.e. throughout the 
time evolution of the fireball), but by hadrons. Due to their strong interac- 
tions, the fireball is opaque to quarks and hadrons for most of its dynamical 
history, and only at the end, after expansion has diluted the matter sufficiently, 
the hadrons decouple. As in cosmology, where, after 400000 years of complete 
opaqueness, the universe became transparent to photons quite suddenly after 
electron-ion recombination, this kinetic decoupling process appears to happen 
in bulk and rather suddenly also in heavy-ion collisions: two-particle Bose- 
Einstein correlation measurements, which give access to the space-time struc- 
ture of the fireball at the point of freeze-outp indicate for Pb+Pb collisions at 
the SPS an emission duration for pions (the most abundantly produced parti- 
cles) of not more than 2-3 fm/c, after a total expansion time of at least 9-10 
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fm/cESLlI Consistently with that, most of the pions_axe emitted from the. bulk 
of the source and not only from a thin surface layerL3'E2l It is also knownE3 that 
in heavy-ion collisions kinetic freeze-out of hadrons is triggered by expansion, 
not by the finite size of the source, again as in the early universe. 

What is, then, the temperature of this thermal hadron radiation? As for 
the cosmic microwave background, it is determined from the energy spectrum 
of the radiated particles and, as in the consmological context, this energy spec- 
trum is affected by (Hubble) expansion. Moreover, as in cosmology, the kinetic 
decoupling temperature reflected in the energy or momentum spectra differs 
from the chemical freeze-out temperature reflected in the particle abundances; 
it is considerably lower. 

Expansion flow affects the observed momentum spectra in two distinct 
ways: since the source flows towards the_d,etector, the single-particle distri- 
butions are flattened by a blueshift effectp^ and as different parts of the ex- 
panding source recede from each other, the homogeneity regions (Hanbury 
Brown/Twiss (HBT) size parameters) measured by Bose-Einstein jmterfero- 
metry (see Sect. |1.2| ) are reduced by velocity gradients in the sourceQcIl 

For the blueshift effect on single-particle transverse mass spectra one must 
distinguish two domains: In the region of relativistic momenta, p±^mQ, the in- 
verse slope Tapp of all particle species is the same and given by the blueshift for- 
mula T'app=Tf W Y^T^iy. This formula does not allow to disentangle the average 

radial flow velocity {vj_) and freeze-out temperature T{. In the non-relativistic 
domain p_L<Cmo the inverse slope is given approximately by T'app=rf-|-TOo(w^), 
and the rest mass dependence of the "appiaxpiit temperature" (inverse slope) 
allows to determine ILand (v^) separatelyE2lE3 (In pp collisions no mo-depen- 
dence of Tapp is seenE3) 

The left diagram in Fig. |6| shows a compilation of measured slope parame- 
ters Tapp for a variety of hadron species in Pb+Pb collisions at the SPS. While 
a rise with the rest mass mo is clearly seen, providing strong evidence for radial 
flow, some of the detailed features lead to ambiguities in the separation of tem- 
perature and flow. The pion slopes are very sensitive to the p^-region in which 
the flts are performed, due to strong resonance decay contributions at low p±. 
Also, pions are alwaysjielativistic such that the formula Tapp=Tf-|-mo(v^) can- 
not be used for thempi and so they don't fit very well into the systematics of 
that Figure. Finally, the il baryons (and perhaps alsolhe S's) exhibit steeper 
slopes than expected from this formula. It was arguecO that this reflects their 
earlier kinetic freeze-out due to an absence of strong scattering resonances with 
the dominating pions; these are essential for the kinetic re-equilibration of the 
other hadron species. 
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A less ambiguous separation of temperature and flow is possible by com- 
bining single-particle m_L-slopes with the pair-momentum dependence of the 
transverse HBT radius extracted from 2-pion correlations (right diagram in 

Fig. H). The latter is approximately given bjo R±{mj_) k. R I Jl + £,{v'\)^^ 
mass dependence of inverse slopes 
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Figure 6: Left: Inverse transverse slope parameter "T" for different hadron species as a fuirri 
tion of their rest masses, measured for Pb+Pb collisions at the SPS by various experimentsEj 
Right: Constraints on possible combinations of kinetic freeze-out temperature T and collec- 
tive radial flow velocity /3^ from the transverse momentum spectra of negative hadrons {h~ ) 
and deuterons (d) and from the transverse momentum dependence of the transverse HBT 
radius parameter Rj_('mj_) (27r-BE)E3 



(see Eq. (H)) where to_i_ is the average transverse mass of the pions in the pair, 
Tf is the kinetic freeze-out temperature, and ^ is a number of order 1 which 
depends on the detailed transverse velocity and density profiles. Since this 
formula and the one for the inverse single-particle slopes provide orthogonal 
correlations between Tf and {vi_) (see Fig. g), their combination leads to a 
rather accurate separation of both. After proper translation of the symbols 
one flndfS for the average radial flow velocity {v±) « 0.5 c and for the tem- 
perature of the thermal hadron radiation Tf « 100 MeV. (The somewhat lower 
Tf-value thanJji Fig. o results from a more accurate fit to the single-particle 
h~ spectrumo) 

Similar analyses were performed for smaller collision systems and at lower 
beam energies; the results are shown by the open symbols in Fig. pi As in 
the Early Universe freeze-out of the momentum spectra is seen to happen later 
than the decoupling of particle abundances, at significantly lower temperatures. 
The ability of the system to remain in a state of approximately local thermal 
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equilibrium while building up collective expansion flow and cooling wa^ recently 
demonstrated in microscopic simulations with the URQMD modelEB 

4-3 Initial energy density 

Two-particle Bose-Einstein correlation measurements give access to both the 
geometry and collective dynamics of the fireball at the point of kinetic freeze- 
outQ'EZl One not only finds signs of strong radial fiow of about 0.5 c, as just 
discussed, but also evidence for strong transverse growth of the fireball be- 
tween impact and freeze-out, by more than a factor 2EJ The two observa- 
tions are dynamically consistent, given the. lower limit on the total expansion 
time Tf >8 fm/c which can be obtaincdEHtJ from the longitudinal HBT ra- 
dius i?|| . Knowing the freeze-out temperature Tf and collective fiow velocity, 
the thermal and collective fiow energy density of the system at freeze-out can 
be calculated. From the measured transverse growth factor and the known 
longitudinal expansion pattern (all extracted from single-particle spectra and 
HBT measurements) one can estimate the total geometric expansion-factor of 
the fireball between the onset of transverse expansion and freczc-outEO Energy 
conservation then gives an estimate of the energy density at the beginning of 
transverse expansion, by multiplying the freeze-out energy density with the 
volume expansion factor. l^ 

This estimate has the advantage over the one using Bjorken's formulalj that 
both factors are determined more or less directly from (single-particle and cor- 
relkation) measurements, thus avoiding uncontrolled model assumptions (like 
the identification of momentum-space with coordinate-space rapidity densities) 
and free parameters (as the equilibration time tq in Bjorken's formula). For 
Pb-|-Pb collisions at the SPS one findsEH in this way initial energy densities 
of 2.5 - 4 GeV/fm'^. This is comfortably above the critical energy density for 
deconfinement, consistent with the other arguments for an early pre-hadronic 
stage given above. Also, this energy density must have been at least partially 
thermalized because pressure (a consequence of thermalization) is necessary to 
drive the transverse expansion. 

4-. 4 Elliptic flow: evidence for ''early pressure" 

In addition to radial transverse flow, which is typical for central collisions, 
two other types of collective flow occur in collisions with finLteJinpact para- 
meter: in-plane directed flow ("bounce-off") and elliptic flowtj'Ej. While the 
first of these is concentrated at forward and backward rapidities, elliptic flow 
is strongest at midrapidity. At very high energies the directed fiow becomes 
weakp3 and only elliptic fiow survives. Both types of directed flow can be 
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identified by a harmonic analysis of tlie azimuthal dependence of the single- 
particle distributions around the beam axis, with iitplane (elliptic) directed 
flow given by the flrst (second) harmonic coefficientLJ 

Elliptic flow arises from the initial elliptic spatial deformation of the over- 
lap region of the two colliding nuclei in the transverse plane. If rescattering be- 
tween secondaries created in this region builds up pressure sufficiently quickly, 
the resulting elliptic anisotropy of the pressure gradients causes an elliptic 
deformation in the developing flow pattern: along the shorter dimension the 
pressure gradient is larger and the flow develops faster. This quickly reduces 
the geometric deformation, i.e. the developing elliptic flow eliminates its own 
cause. This is the main reason why elliptic flow is particularly sensitive to the 
pressure and equation of state in the early stages of the collijaion whereas radial 
flow receives strong contributions also from the late stagesB 




3 4 

rapidity 



Figure 7: Second harmonic coefficient V2 of the azimuthal distribution (elljptic flow) of 
protons in Pb+Pb collisions at the SPS, as a fpjction of rapidity. NA49 datall3 are shown 
together with RQMD and VENUS simulations^ 

Elliptic fkov was measured in Pb+Pb collisions at the SPS by the NA49 
collaborationEJ Figure shows the data -for protons (due to their larger mass 
they show a stronger signal than pionscil), together with simulations using 
RQMD and VENUS event generators. Clearly the latter underpredict the effect 
significantly. This failure is due to the particular way these codes parametrize 
particle production: colliding nucleons create strings which, after a formation 
time of about 1 fm/c, decay directly into hadrons. During the initial stage 
of highest energy density the energy is thus stored in non-interacting strings 
which do not contribute to the pressure. In other words, their equation of state 
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is "ultra-soft" la 

Sorge has recently shownQ that a modification of the RQMD code which 
simulates a harder QGP equation of state during the early stage is able to 
generate the elliptic flow signal. Moreover, he found that the entire effect is 
created during the high density stage with e > ecr;|-5tfter hadronization the 
elliptic flow quickly saturatesil The elliptic flow datacd can thus be viewed as 
a rather direct glimpse of the quark-gluon plasma. 

4-5 Absence of non- statistical event-by-event fluctuations 

In the introduction I mentioned that in cosmology the only direct signature 
from the time before chemical and thermal freeze-out seem to be the recently 
discovered small temperature fluctuations in the microwave background. Gen- 
erated by gravitation, they have no analogue in heavy-ion physics. There, 

NA49 Pb+Pb Event-by-Event Fluctuations 
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Figure 8: Event-by-event fluctuations of (p±) and K/n in Pb-|-Pb collisions at the SPStJ 
The solid lines denote the statistical fluctuations from mixed events. 

however, we can study large numbers of collision events and investigate the 
fluctuations of physical observables from event to event. Such fluctuations are 
expected for statistical reasons because of the flnite size of the collision system 
and the finite number of produced particles. If the collision fireballs are truely 
thermalized, finite number effects and quantum statistical fluctuations should 
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be the only reason for event-by-event fluctuationsEJ 

Figure shows the measured distributions for the average transverse mo- 
nicntiini {p±) and for the K/ir ratio in individual Pb-|-Pb coUision events at the 
SPSe3 These distributions are perfect Gaussians to a level of a few times 10~^; 
moreover, the widths of these Gaussians agree very accurately with expecta- 
tions from mixed events, i.e. the fluctuations in the measured quantities from 
event to event are consistent with purely statistical (e.g. thermal) fluctuations. 
No signs for large "critical fluctuations" expected near a phase transition or 
of other dynamical fluctuations are visible. If, as argued above, the chemical 
composition, which leads to a given average K/tt ratio, and the thermaliza- 
tion and flow pattern, which causes a certain average transverse momentum of 
the emitted hadrons, reflect interesting physics of the fireball evolution, every 
Pb+Pb event at the SPS reflects it in exactly the same way, at the level of less 
than a permille! 

5 Conclusions 

I have presented arguments that in heavy-ion collisions at the SPS we have 
seen THE LITTLE BANG: the measured soft hadron spectra show con- 
vincing evidence for Hubble-like (3-dimcnsional) flow, for thermal (hadron) 
radiation, and for "primordial" hadrosynthesis out of a pre-hadronic state 
without quark confinement. There are strong signs for the existence of a pre- 
hadronic stage with non-trivial dynamics, resulting in a global strangeness 
enhancement by about a factor 2, in a particularly strong enhancement of 
multistrange (anti)baryons consistent with a statistical hadronization picture, 
in the suppression of J/ip, Xc a-nd tp' states, and in elliptic fiow. The absence 
of non-statistical fluctuations from event to event indicates that all collision 
events are similar, that these interesting phenomena occur in every event, and 
that decoupling did not take place close to a phase-transition. For chemical 
decoupling (which I argued to happen directly after hadron formation) this im- 
plies that the hadronization process cannot be viewed as an equilibrium phase 
transition; thermal freeze-out was found to occur at much lower temperatures, 
in safe distance from any phase transition. 

The abundance and quality of hadron data now available allows for a 
detailed and quantitative description of the final state in heavy-ion collisions. 
This begins to severely constrain dynamical extrapolations backward in time. 
Although no direct signatures from the pre-hadronic stage have so far been 
seen, the soft hadron data leave little room for doubt that it is there. Whether 
we can call it a quark-gluon plasma which satisfies the criteria spelled out in 



Sec. 1.2 will only become clear at higher collision energies where such direct 
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signals become accessible. At the moment we only know that the pre-hadronic 
state exerts pressure, so a certain degree of thermalization among the partons 
must have occurred. The other important open question is, of course: where 
is the threshold for deconfinement and how do we identify it experimentally? 
Clearly, much work remains to be done but now we know that we are on the 
right track. 
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